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The search for topological superconductors which support Majorana fermion exci¬ 
tations has been an important topic in condensed matter physics. In this work, we 
propose a new experimental scheme for engineering topological superconductors. In 
this scheme, by manipulating the superlattice structure of organic molecules placed on 
top of a superconductor with Rashba spin-orbit coupling, topological superconducting 
phases can be achieved without fine-tuning the chemical potential. Moreover, super¬ 
conductors with different Chern numbers can be obtained by changing the superlattice 
structure of the organic molecules. 


I. INTRODUCTION 


Topological superconductors are new states of mat¬ 
ter which are gapped in the bulk and support gapless 
Majorana fermion excitations at the boundary HHS]. It 
has been pointed out that Majorana fermions are self- 
Hermitian particles and obey non-Abelian statistics with 
potential applications in fault-tolerant quantum compu¬ 
tations IMI] and in spintronics PIT3] . Due to the ex¬ 
otic properties of Majorana fermions, the search for topo¬ 
logical superconductors has been an important topic in 
recent years. 

It was first pointed out that Majorana fermions can ex¬ 
ist in chiral p-wave superconductors where spin-polarized 
electrons are paired to form Cooper pairs Enn. It was 
shown by Kitaev that the essence of creating chiral p- 
wave superconductor is to pair electrons from an odd 
number of partially occupied subbands at the Fermi en¬ 
ergy 0- However, intrinsic chiral p-wave superconduc¬ 
tors are yet to be identified. 

After the discovery of topological insulators, Fu and 
Kane proposed that conventional s-wave superconduc¬ 
tors in proximity to a topological insulator surface can 
be used to engineer p-wave topological superconductors 
m- Their proposal was based on the fact that: 1) an 
odd number of partially occupied surface bands exist on 
the surface of topological insulators, and 2) the surface 
electrons are helical, in the sense that electrons with op¬ 
posite momentum have opposite spin due to spin orbit 
couplings (SOC). As a result, superconducting pairings 
can be effectively induced by an s-wave superconductor 
through proximity effect. The resulting s-wave supercon¬ 
ductor can be mapped to an effective p-wave topological 
superconductor. This proposal has been pursued by sev¬ 
eral experimental groups [TBHTH] but it is a challenge to 
distinguish the Majorana fermions from other low energy 
fermionic excitations in the system. 

Another promising way to create a chiral p-wave su¬ 
perconductor is to apply a Zeeman field to a semiconduc¬ 
tor with Rashba spin-orbit coupling (SOC) and induce 
s-wave pairing through proximity effects [TqHZS]. First, 
the Zeeman field opens an energy gap at the F point such 



FIG. 1: Schematic ID band structures (a) In the presence of 
Rashba SOC and a Zeeman field. When the chemical poten¬ 
tial jjL is in the Zeeman gap, the outer subband is partially 
filled, (b) In the presence of a Zeeman field, small Rashba 
SOC and an infinitesimal superlattice potential. The super¬ 
lattice potential with period n/kp reduces the Brillouin zone 
to [—kF^kp]- (c) With finite super lattice potential, an energy 
gap opens at the Brillouin zone boundary and results in a 
single partially occupied band at the Fermi energy. Induc¬ 
ing superconductivity at the Fermi energy will result in an 
effective p-wave superconductor. 


that there is only one partially filled subband at Fermi 
energy if the Fermi energy is inside the Zeeman gap as 
depicted in Fig. la. Second, the partially occupied sub¬ 
bands have helical properties due to Rashba SOC. As a 
result, pairing of the helical electrons can also be induced 
by an s-wave superconductor and result in an effective 
spinless p-wave superconductor. Experiments following 
this scheme have been performed but the results are not 
yet conclusive [MEU- One of the most challenging parts 
of the experiment is to fine tune the chemical potential 
such that it falls within the Zeeman gap at the F point, 
where the size of the Zeeman gap is in the order of meV. 

In this work, we propose that putting superlattices on 
top of a system with Rashba split bands can be used 
to realise effective p-wave superconductors without fine- 
tuning the chemical potential. To illustrate the ideas, we 
first consider a one dimensional system with Rashba split 
bands and where the Fermi energy is far away from the 
band bottom as illustrated in Fig.lb. We demonstrate 
how a superlattice potential can open an energy gap for 
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one of the spin subbands at the Brillouin zone boundary 
such that there is only one partially occupied helical sub¬ 
band at the Fermi energy, as depicted in Fig.lc. Inducing 
s-wave superconductivity on the system will result in a 
topological superconductor. 

Second, we consider a two dimensional system with 
Rashba split bands. We show that superlattices, together 
with time-reversal symmetry breaking, change the Chern 
number of certain subbands of the system and make the 
system topological. Interestingly, we show that states 
with different Chern numbers can be engineered using 
different superlattice structures and different lattice spac¬ 
ing. Square, rectangular and triangular superlattices are 
studied in detail as depicted in Fig.3 and Fig.4. Im¬ 
portantly, recent experiments have demonstrated that 
organic molecules can induce superlattice potentials on 
metal surfaces to create different band structures such as 
Dirac band structures IMES]. We believe that inducing 
superlattice potential and s-wave superconductor on An 
or Bi surfaces which possess Rashba surface bands, can 
be used to create topological superconductors without 
fine-tuning the chemical potential. 

II. ID REALIZATION 

Schematic illustration — To illustrate the importance 
of band folding introduced by superlattices, we first con¬ 
sider a wire with Rashba spin-orbit coupling. In the basis 
of the Hamiltonian of the wire in the presence 

of a magnetic field and s-wave pairing can be written as 

[I9H23]: 



FIG. 2: (a) The schematic illustration of a ID wire with the 
potential of every third site is modified, (b) The spectral 
function of the last site of Hidtb with A = 0, t = 1, an = 
O.Olt and Vx = 0.05t. (c) The spectral function A(ct;) with 
parameters as in (b) except A = 0.015t. (d) The normal 

state band structure with an = 0.2t and Vx = 0.05t. (e) 

The topological invariant Nbdi of Hidtb- The topological 
regimes in (c) and (e) coincide. 


^2 

Hid = + VxTz(^x + OLRkTzcry + Aryay ( 1 ) 

Here, m is the effective mass of the electrons, aR is the 
Rashba SOC strength, Vx is the Zeeman energy from an 
external magnetic field along the wire. As pointed out 
previously, the wire becomes a topological superconduc¬ 
tor when the chemical potential is inside the Zeeman gap 
such that |/i| < — A^, as depicted in Fig. la. In 

this case, the chemical potential intercepts a single non¬ 
degenerate band and the system can be mapped to a 
spinless p-wave superconductor [ToHlS]. However, it is 
very difficult to gate the chemical potential to the Zee- 
man gap experimentally. 

In this section, we point out that the system can be¬ 
come a topological superconductor without fine-tuning 
the chemical potential. The key is to introduce a super¬ 
lattice potential u{x) with period approximately equal to 
211 jkR where the kp is the Fermi wave vector. 

To be specific, we first discuss the case with weak 
Rashba SOC such that aRkp <C 14 <C t. These parame¬ 
ters are applicable to realistic semiconducting wires used 
in recent experiments [HHHEn]. In this case, we have 
two subbands separated by the Zeeman energy. When 
the superlattice potential is infinitesimal, the size of the 
Brillouin zone is reduced such that T/ci? become the 


zone boundary points as illustrated in Fig.lb. However, 
when the superlattice potential is finite, it couples states 
with the same spin at k and /c ± G, where G = 2k f 
and ug is the corresponding Fourier component of u{x). 
This can be easily seen in the second quantized form 
whereby the superlattice potential term can be written 
as UG^i^1pk±Ga- 

This coupling will result in an energy gap at the re¬ 
duced Brillouin zone boundary as depicted in Fig.lc. 
Denote I'lpik)) as an eigenstate of HiD{k) with A = 
0, the energy gap at the reduced zone boundary 

IS (' 0 ( f )\ g|' 0 ( f )) W+q,2^0274+^0, 

From Fig.lc, we see that if the Fermi energy is within 
certain energy gaps of the Brillouin zone boundary, there 
is only one partially filled subband at the Fermi energy. 
Further introducing an s-wave pairing into the system 
will result in a topological superconductor. It is impor¬ 
tant to note that the key to reach the topological regime 
is to introduce a superlattice potential with appropriate 
lattice constant. 


Tight-binding calculations — To verify the above 
claims, we construct a tight-binding model to describe 
Hid and the superlattice potential. To be specific, the 
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superlattice potential is simulated by changing the on¬ 
site energy of one site out of every three sites of the 
original wire as illustrated in Fig. 2a. The model can be 
written as: 

HidTB = Y.- + h.C.) - II-lpl 

+ Y- + «4i,a^3i,a +/l-C. (2) 

^ —' 2 , 0 ; ’ ’ ’ 

Here, i denotes the lattice sites and the spin indices are 
denoted by a and /3, t and fi are the hopping amplitude 
and the chemical potential respectively, aR is the Rashba 
coupling strength, is the Zeeman energy, A is the s- 
wave pairing amplitude, and u is the on-site potential 
which is non-zero for every three sites. 

As a result of band folding, the new Brillouin zone 
boundary is located at ±7r/3a where a is the lattice spac¬ 
ing of the tight-binding model. The energy spectrum of 
the model in the absence of superconducting pairing is 
shown in Fig.2b. As expected, energy gaps open at the 
Brillouin zone boundary as well as the F point. The sys¬ 
tem is expected to be a topological superconductor when 
the chemical potential is within the regime where an odd 
number of subbands are partially occupied. 

To verify the existence of the Major ana end states 
of HiDTBi we plot the spectral function A(co’) = 
Tr[/mG(R, R, cj)] of the last site R of a semi-infinite 
wire as a function of chemical potential and energy in 
Fig.2c. Here, G is the lattice Green function of Hidtb 
1ST]. From the spectral function, it is evident that zero 
energy modes, which correspond to Majorana end states 
of the wire, appear at the regime where an odd number of 
subbands are partially occupied. The topological regime 
at low chemical potential is due to the opening of the 
Zeeman gap near the band bottom which would appear 
even in the absence of the superlattice potential. On the 
other hand, the topological regime created by the open¬ 
ing of an energy gap at the Brillouin zone boundary at 
the Fermi energy is due to the superlattice potential. 

Topological invariant — To further study the topolog¬ 
ical properties of the system, we write Hidtb in the 
momentum basis. Due to the superlattice potential, 
there are three sublattices in each unit cell. The ba¬ 
sis can be written as where 

^/c,a = ['0/c,a(l),'0/c,a(2),'0/c,a(3)] and 'ipk,a{j) dcnotcs an 
electron annihilation operator of the j sublattice. In this 
basis, Hidtb is in the BDI class due to the presence 
of the particle-hole symmetry P = r^K and the chiral 
symmetry C = where K is the complex conjugate 
operator and operates on the particle-hole sector. As 
a result, we can classify Hidtb by a topological invariant 
Hbdi where Nbdi is defined in the Method section. The 
topological invariant Nb di as a function of the chemical 
potential is shown in Fig.2e. By comparing the topo¬ 
logical regime in Fig.2e with the gap opening regime at 


Brillouin zone boundary in Fig.2b, it is evident that the 
system is topologically non-trivial when an odd number 
of subbands are partially occupied. 

We note that introducing a superlattice potential can 
also create topological regimes at the Fermi energy when 
Rashba energy is strong. For example, the normal state 
spectrum for Hidtb with strong Rashba SOC such that 
OiRkr ^ 14 is illustrated in Fig.2d. In this case, due to 
the large Fermi momentum difference between the two 
subbands, the lower subband with larger is gapped 
out at the Brillouin zone boundary by the superlattice 
potential. On the other hand, the upper subband remains 
gapless. The system is topologically non-trivial when an 
odd number of subbands in the normal state are partially 
occupied and when s-wave superconductivity is induced. 


III. SUPERLATTICE IN 2D SYSTEMS 

The above section illustrated how ID superlattices 
with suitable lattice spacing can be used to engineer a ID 
topological superconductor. In this section, we demon¬ 
strate that 2D superlattices on metal surface states can 
be used to engineer topological superconductors with dif¬ 
ferent Chern numbers without fine-tuning the chemical 
potential. 

Our proposal is inspired by a recent experiment in 
which coronene molecules were deposited on top of cop¬ 
per [111] surfaces [29]. The superlattice potential induced 
by the molecules creates a honeycomb lattice structure 
on the surface and results in Dirac spectrum of the sur¬ 
face states. The positions of the coronene molecules on 
the copper surface can be manipulated by the STM tips 
so that the lattice spacing and structures of the superlat¬ 
tice potential are highly tuneable. Several types of arti¬ 
ficially made defect states and zero energy states associ¬ 
ated with the zig-zag edge of the sample, were observed, 
demonstrating the great tuneability of the setup. With 
this recent advancement in experimental technique, we 
believe that desirable superlattice potential can be en¬ 
gineered on metal surfaces with strong Rashba surface 
states as well. 

Square lattice — To be specific, we consider Au(lll) 
surfaces and the Au thin layer is placed on top of a su¬ 
perconducting substrate as illustrated in Fig.3a. It is 
shown both theoretically and through ARPES experi¬ 
ments that the surface states have large Rashba splitting 
of about lOOmeV [32]. It was proposed that Au(lll) 
surface states, when combined with Zeeman field and su¬ 
perconductivity, can be used to generate topological su¬ 
perconductors. However, the chemical potential has to 
be tuned to the Zeeman gap near the Rashba band bot¬ 
tom [33]. Since the Fermi energy of Au is about 0.4eV 
above the Rashba band bottom [32], it is very difficult 
to gate the chemical potential to the Zeeman gap. As 
we show below, a suitable superlattice potential and a 
Zeeman field can make the Au thin film a topological 
superconductor. 
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FIG. 3: (a) Organic molecules on the top of An (111) sur¬ 

face where An is grown on top of an s-wave superconductor, 
(b) The outer (green) and inner Fermi (blue) circles in the 
superlattice Brillouin zone. The square indicates in reduced 
Brillouin zone due to superlattice potential, (c) The band 
closing and reopening at X point, as a function of 14. 14c 
denotes the critical value of 14 at which the band gap closes, 
(d) The energy spectrum of An surface with square super¬ 
lattice potential. The tight-binding model H 2 dtb used is 
described in the Method section. The superlattice is chosen 
such that the Chern number is -2. It is evident that there are 
two crossings at zero energy corresponding to the two zero 
energy Major ana modes at each edge, (e) The Chern num¬ 
ber of H 2 DTB as a function of chemical potential for a fixed 
superlattice configuration. 


We first place organic molecules on top of the Au(lll) 
surface to induce a superlattice potential as illustrated in 
Fig.3a. If the superlattice potential has a square struc¬ 
ture with lattice spacing Tr/Zci?, where kp is the Fermi 
wave vector of the Rashba band with larger Fermi mo¬ 
mentum, the X and Y points (T/ci?, 0), (0, T/c^) becomes 
the Brillouin zone boundary points as depicted in Fig.3b. 
From ARPES experiments, kp is about 0.2/A [32] and 
the superlattice spacing should be about 1.6nm. This 
superlattice spacing falls within the experimentally ac¬ 
cessible regime. 

Now, if an s-wave pairing is induced on Au(lll) sur¬ 
face, a pairing gap opens at the Fermi energy as illus¬ 
trated in Fig.3c. The energy spectrum in Fig.3c is based 
on a tight-binding model introduced in the Method sec¬ 
tion. Due to the pairing gap, one may classify this two 
dimensional system by Chern numbers. However, due to 
time-reversal symmetry, the total Chern number of the 
occupied bands have to be zero. Indeed, from the tight- 



kb/jT |Li/t 


FIG. 4: (a) The schematic illustration of organic molecules, 

denoted by the disks, on An (111) surface. The local po¬ 
tential energy of the Au under the molecules are modified. A 
hexagonal unit cell with 16 Au sites, used in the tight-binding 
model, is shown, (b) The hexagon denotes the reduced Bril¬ 
louin zone due to the superlattice potential. The superlattice 
is chosen such that the M points touch the Fermi surface of 
the outer Rashba band, (c) The energy spectrum of a su¬ 
perconductor with Ghern number 3. The superconductor is 
obtained by adding superconductivity and a Zeeman field to 
a system with reduced Brillouin zone given in (b). (d) The 
Ghern number as a function of chemical potential for a fixed 
superlattice configuration. 


binding model calculations, one can show that the outer 
and inner occupied bands have Chern number 1 and -1 
respectively. 

Importantly, a Zeeman field can close the bulk gap at 
the time-reversal invariant points X and Y and reopen it 
again as illustrated in Fig.3c. After the gap reopening, 
the Chern number of the outer subband is changed by 
two, from 1 to -1. As a result, the total Chern number of 
the occupied bands is -2. The energy spectrum of the sys¬ 
tem with open boundary condition in the y-direction and 
periodic boundary condition in the x-direction is shown 
in Fig.3e. It is evident that in-gap Majorana edge states 
appear in the energy spectrum. The tight-binding model 
used is described in the Method section. 

Interestingly, if the superlattice potential is rectangu¬ 
lar such that only the energy gap at X point in Fig.3b 
is closed and reopened by the Zeeman field, the Chern 
number of the outer band is zero. As a result, the total 
Chern number of the system is -1 instead. 

Triangular lattice — One of the great advantages of the 
superlattice scheme is that many different lattice struc¬ 
tures can be constructed on the metal surfaces. An in¬ 
teresting possibility is to construct a triangular lattice on 
top of Au(lll) surface as depicted in Fig.4a. If the super¬ 
lattice spacing is approximately it/ kp^ the new Brillouin 






















































5 


zone will have hexagonal shape and the three M points 
will touch the outer Rashba bands of the Au(lll) sur¬ 
face states as shown in Fig.4b. Adding an s-wave pairing 
potential will gap out the Fermi surfaces and the total 
Chern number of the occupied bands is zero. However, 
increasing the Zeeman field can close and re-open the gap 
at the three M points such that the Chern number of the 
outer band is changed by 3 from 1 to -2. As a result, the 
total Chern number of the occupied bands is -3. This 
results in three chiral edge states at a zig-zag edge of 
the sample as shown in Fig.4c. The Chern number as 
a function of chemical potential for a fixed superlattice 
configuration is shown in Fig.4d. 


IV. DISCUSSION 

In the above sections, we discussed how non-magnetic 
superlattice potentials and uniform magnetic field can be 
used to engineer topological superconductors in systems 
with Rashba split bands. In particular, we pointed out 
that the superlattice potential on Au(lll) surface can 
be induced by non-magnetic coronene molecules. How¬ 
ever, a strong external uniform magnetic field is not nec¬ 
essary. Indeed, the non-magnetic molecules can be re¬ 
placed by organic molecules with magnetic centers. For 
example, it has been shown that Co-phthalocyanine and 
Fe-porphyrin adsorbed on Au(lll) surface carries mag¬ 
netic moments and exhibit Kondo resonances [34ll36] . As 
a result, the moments can be aligned by a weak external 
magnetic field. These ferromagnetic superlattice poten¬ 
tial can result in topological superconductors similar to 
uniform magnetic field with non-magnetic superlattices. 

It is also important to note that the superlattice 
method proposed in this work is very general. The 
Au(lll) surface can be replaced by other 2D systems 
such as Bi (111) [37] and Sb (111) |38l39j surfaces which 
also support Rashba bands. 

Another possible way to realize topological supercon¬ 
ductors is to place magnetic molecules on superconduct¬ 
ing Pb thin films. It has been shown that Pb thin films 
on Si substrate possess Rashba split bands m and be¬ 
come superconducting at low temperatures |41l42j . Im¬ 
portantly, it has been shown that magnetic metalorganic 
molecules on Pb surfaces can induce Yu-Shiba states m- 
Therefore, we believe that superlattice potentials on Pb 
surfaces induced by magnetic molecules can be used to 
realize topological superconductors without fine-tuning 
the chemical potential as demonstrated above. 


V. METHOD 


With the matrix A{k), one can define the phase 6>, 

z{k) = = Det[A(fc)]/|Det[A(fc)]|, (4) 

such that 0{k) = utt at /c = 0, Ttt with integer n. The 
winding number of 0{k) can be used as the topological 
invariant which characterizes the Hamiltonian HiDTB{k) 
|44l45j . The winding number NbdIi which counts the 
number of Majorana end states at one end of the wire, 
can be written as 


dz{k) 


(5) 


Nbdi as a function of chemical potential is shown in 
Fig.2e. When comparing Fig.2b to Fig.2e, it is evident 
that the system is topological when the chemical poten¬ 
tial intercepts a single subband. 

2D tight-binding models - The Rashba bands of the 
Au(lll) surface with the superlattice potential is de¬ 
scribed by a tight-binding model. The lattice structure 
of An (111) surface is triangular as depicted in Fig.4a. To 
include the effect of the superlattice, we choose a hexago¬ 
nal unit cell with 16 sites and the on-site chemical poten¬ 
tial of the 7 sites inside the muffin-tin circles are modified. 
The tight-binding model can be written as: 


H 2 DTB = Y M Y +«(*) Y dV’i 

{i,j) i i 

- iaR Y da ^ 

- v Y (dtV’it - ^4^4) - ^ Y dtd + h.c. 

i i 

(6) 

Here, is a unit vector connecting sites i and j, Bz is 
a unit vector in the 2 : direction, u (i) is nonzero for sites 
inside the muffin-tin circle as depicted in Fig.4a. 

In the square lattice case, for simplicity, we have cho¬ 
sen a square lattice tight-binding model to simulate the 
Rashba bands. A unit cell of 16 sites is used and the su¬ 
perlattice is simulated by changing the local chemical po¬ 
tential of some of the sites in the unit cell, similar to the 
triangular case. The square lattice model can be justified 
by the fact that the Fermi surfaces of the Rashba bands 
of An (111) surface are circular and far from the original 
Brillouin zone boundary. The same circular Fermi sur¬ 
face can be reproduced by square lattice models as long 
as the Fermi surface is far away from the Brillouin zone 
boundary. 


Calculation of Nbdi — To calculate the ID topologi¬ 
cal invariant of HidtBi we note that in the basis which 
diagonalises (7, Hidtb can be block diagonalised as: 

H^nrnik) = ( ^ (3) 
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the finishing stage of this work, we noted that there are portance of band folding was not emphasised in these 
proposals using 2D arrays of Yu-Shiba states to engineer works, 
topological superconductors even though the im- 
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